(* Content-type: application/mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 7.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       145,          7]
NotebookDataLength[     30198,        875]
NotebookOptionsPosition[     28504,        825]
NotebookOutlinePosition[     28899,        842]
CellTagsIndexPosition[     28856,        839]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[TextData[StyleBox["\:0424\:0443\:043d\:043a\:0446\:0438\:0438 \:043f\
\:043e\:0441\:0442\:0440\:043e\:0435\:043d\:0438\:044f \:041b\:0410\:0427\
\:0425 \:0438 \:043f\:0435\:0440\:0435\:0445\:043e\:0434\:043d\:043e\:0439 \
\:0444\:0443\:043d\:043a\:0446\:0438\:0438",
 FontFamily->"Times New Roman",
 FontSize->48,
 FontColor->RGBColor[0., 0., 1.]]], "Text"],

Cell[BoxData[
 RowBox[{"Clear", "[", "\"\<Global`*\>\"", "]"}]], "Input"],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"Bode", "[", 
    RowBox[{"\[CapitalPhi]_", ",", "f0_", ",", "fn_"}], "]"}], ":=", 
   RowBox[{"Module", "[", 
    RowBox[{
     RowBox[{"{", "\[Omega]", "}"}], ",", "\[IndentingNewLine]", 
     RowBox[{"GraphicsColumn", "[", "\[IndentingNewLine]", 
      RowBox[{"{", "\[IndentingNewLine]", 
       RowBox[{
        RowBox[{"LogLinearPlot", "[", 
         RowBox[{
          RowBox[{"20", "*", 
           RowBox[{"Log10", "[", 
            RowBox[{"Abs", "[", 
             RowBox[{
              RowBox[{"(", 
               RowBox[{"\[CapitalPhi]", "/", "p"}], ")"}], "/.", 
              RowBox[{"{", 
               RowBox[{"p", "\[Rule]", 
                RowBox[{"I", "*", "\[Omega]"}]}], "}"}]}], "]"}], "]"}]}], 
          ",", 
          RowBox[{"{", 
           RowBox[{"\[Omega]", ",", 
            RowBox[{"10", "^", "f0"}], ",", 
            RowBox[{"10", "^", "fn"}]}], "}"}], ",", 
          RowBox[{"AxesOrigin", "\[Rule]", 
           RowBox[{"{", 
            RowBox[{
             RowBox[{"10", "^", "f0"}], ",", "0"}], "}"}]}]}], "]"}], ",", 
        "\[IndentingNewLine]", 
        RowBox[{"LogLinearPlot", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"180", "/", "Pi"}], "*", 
           RowBox[{"Arg", "[", 
            RowBox[{
             RowBox[{"(", 
              RowBox[{"\[CapitalPhi]", "/", "p"}], ")"}], "/.", 
             RowBox[{"{", 
              RowBox[{"p", "\[Rule]", 
               RowBox[{"I", "*", "\[Omega]"}]}], "}"}]}], "]"}]}], ",", 
          RowBox[{"{", 
           RowBox[{"\[Omega]", ",", 
            RowBox[{"10", "^", "f0"}], ",", 
            RowBox[{"10", "^", "fn"}]}], "}"}], ",", 
          RowBox[{"AxesOrigin", "\[Rule]", 
           RowBox[{"{", 
            RowBox[{
             RowBox[{"10", "^", "f0"}], ",", "0"}], "}"}]}]}], "]"}], ","}], 
       "}"}], "\[IndentingNewLine]", "]"}]}], "]"}]}], 
  "\[IndentingNewLine]"}]], "Input"],

Cell[BoxData[
 RowBox[{
  RowBox[{"Step", "[", 
   RowBox[{"\[CapitalPhi]_", ",", "tn_"}], "]"}], ":=", 
  RowBox[{"Module", "[", 
   RowBox[{
    RowBox[{"{", "t", "}"}], ",", 
    RowBox[{"Plot", "[", 
     RowBox[{
      RowBox[{"InverseLaplaceTransform", "[", 
       RowBox[{
        RowBox[{"\[CapitalPhi]", "/", "p"}], ",", "p", ",", "t"}], "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"t", ",", "0", ",", "tn"}], "}"}], ",", 
      RowBox[{"PlotRange", "\[Rule]", "Full"}]}], "]"}]}], "]"}]}]], "Input"],

Cell[TextData[StyleBox["\:041f\:0440\:043e\:0434\:043e\:043b\:044c\:043d\:043e\
\:0435 \:0434\:0432\:0438\:0436\:0435\:043d\:0438\:0435",
 FontFamily->"Times New Roman",
 FontSize->48,
 FontColor->RGBColor[0., 0., 1.]]], "Text"],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
   "\:0412\:0435\:043a\:0442\:043e\:0440", " ", "X", " ", 
    "\:043f\:0440\:043e\:0434\:043e\:043b\:044c\:043d\:043e\:0433\:043e", " ",
     "\:0434\:0432\:0438\:0436\:0435\:043d\:0438\:044f"}], "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{"Xprod", "=", 
   RowBox[{"{", " ", 
    RowBox[{
     RowBox[{"\[Omega]z", "[", "t", "]"}], ",", " ", 
     RowBox[{"\[Alpha]", "[", "t", "]"}], ",", " ", 
     RowBox[{"H", "[", "t", "]"}], ",", " ", 
     RowBox[{"\[CurlyTheta]", "[", "t", "]"}]}], "}"}]}]}]], "Input"],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
   "\:041c\:0430\:0442\:0440\:0438\:0446\:0430", " ", "\:0410", " ", 
    "\:043f\:0440\:043e\:0434\:043e\:043b\:044c\:043d\:043e\:0433\:043e", " ",
     "\:0434\:0432\:0438\:0436\:0435\:043d\:0438\:044f"}], "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"A\[Omega]z", "=", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"c1", "+", "p"}], ",", 
       RowBox[{"-", "1"}], ",", "0", ",", 
       RowBox[{"-", "1"}]}], "}"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"A\[Alpha]", "=", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{
        RowBox[{"c5", " ", "p"}], "+", "c2"}], ",", " ", 
       RowBox[{"p", "+", "c4"}], ",", "c6", ",", "0"}], "}"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"AH", "=", 
     RowBox[{"{", 
      RowBox[{"0", ",", "0", ",", "p", ",", "0"}], "}"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"A\[CurlyTheta]", "=", 
     RowBox[{"{", 
      RowBox[{"0", ",", "0", ",", 
       RowBox[{"-", "c6"}], ",", "p"}], "}"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{"MatrixForm", "[", 
    RowBox[{"Aprod", "=", 
     RowBox[{"Transpose", "[", 
      RowBox[{"{", " ", 
       RowBox[{
       "A\[Omega]z", ",", " ", "A\[Alpha]", ",", " ", "AH", ",", 
        "A\[CurlyTheta]"}], "}"}], "]"}]}], "]"}]}]}]], "Input"],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
   "\:0412\:0435\:043a\:0442\:043e\:0440", " ", "U", " ", 
    "\:043f\:0440\:043e\:0434\:043e\:043b\:044c\:043d\:043e\:0433\:043e", " ",
     "\:0434\:0432\:0438\:0436\:0435\:043d\:0438\:044f"}], "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{"Uprod", "=", 
   RowBox[{"{", 
    RowBox[{"\[Delta]\:0432", "[", "t", "]"}], "}"}]}]}]], "Input"],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
   "\:041c\:0430\:0442\:0440\:0438\:0446\:0430", " ", "B", " ", 
    "\:043f\:0440\:043e\:0434\:043e\:043b\:044c\:043d\:043e\:0433\:043e", " ",
     "\:0434\:0432\:0438\:0436\:0435\:043d\:0438\:044f"}], "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{"MatrBprod", "=", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "c3"}], ",", "0", ",", "0", ",", "0", ",", "0"}], 
    "}"}]}]}]], "Input"],

Cell[BoxData[
 RowBox[{"\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{
   "\:0423\:0440\:0430\:0432\:043d\:0435\:043d\:0438\:044f", " ", 
    "\:043f\:0440\:043e\:0434\:043e\:043b\:044c\:043d\:043e\:0433\:043e", " ",
     "\:0434\:0432\:0438\:0436\:0435\:043d\:0438\:044f"}], "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{"MatrixForm", "[", 
   RowBox[{"prod", "=", 
    RowBox[{"Table", "[", " ", 
     RowBox[{
      RowBox[{
       RowBox[{"Sum", "[", 
        RowBox[{
         RowBox[{
          RowBox[{"Aprod", "[", 
           RowBox[{"[", 
            RowBox[{"i", ",", "j"}], "]"}], "]"}], "*", 
          RowBox[{"Xprod", "[", 
           RowBox[{"[", "j", "]"}], "]"}]}], ",", 
         RowBox[{"{", 
          RowBox[{"j", ",", "1", ",", "4"}], "}"}]}], "]"}], "\[Equal]", 
       RowBox[{
        RowBox[{"MatrBprod", "[", 
         RowBox[{"[", "i", "]"}], "]"}], "*", 
        RowBox[{"Uprod", "[", 
         RowBox[{"[", "1", "]"}], "]"}]}]}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", "1", ",", "4"}], "}"}]}], " ", "]"}]}], 
   "]"}]}]], "Input"],

Cell[TextData[StyleBox["\:0411\:043e\:043a\:043e\:0432\:043e\:0435 \:0434\
\:0432\:0438\:0436\:0435\:043d\:0438\:0435",
 FontFamily->"Times New Roman",
 FontSize->48,
 FontColor->RGBColor[0., 0., 1.]]], "Text"],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
   "\:0412\:0435\:043a\:0442\:043e\:0440", " ", "X", " ", 
    "\:0431\:043e\:043a\:043e\:0432\:043e\:0433\:043e", " ", 
    "\:0434\:0432\:0438\:0436\:0435\:043d\:0438\:044f"}], "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{"Xbok", "=", 
   RowBox[{"{", " ", 
    RowBox[{
     RowBox[{"\[Omega]x", "[", "t", "]"}], ",", 
     RowBox[{"\[Omega]y", "[", "t", "]"}], ",", 
     RowBox[{"\[Beta]", "[", "t", "]"}], ",", 
     RowBox[{"\[Gamma]", "[", "t", "]"}], ",", 
     RowBox[{"\[Psi]", "[", "t", "]"}]}], "}"}]}]}]], "Input"],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
   "\:041c\:0430\:0442\:0440\:0438\:0446\:0430", " ", "\:0410", " ", 
    "\:0431\:043e\:043a\:043e\:0432\:043e\:0433\:043e", " ", 
    "\:0434\:0432\:0438\:0436\:0435\:043d\:0438\:044f"}], "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"A\[Omega]x", "=", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"p", "+", "b1"}], ",", "b6", ",", 
       RowBox[{"-", "b7"}], ",", 
       RowBox[{"-", "1"}], ",", "0"}], "}"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"A\[Omega]y", "=", 
     RowBox[{"{", 
      RowBox[{"a6", ",", 
       RowBox[{"p", "+", "a1"}], ",", 
       RowBox[{"-", "1"}], ",", "0", ",", 
       RowBox[{"-", "1"}]}], "}"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"A\[Beta]", "=", 
     RowBox[{"{", 
      RowBox[{"b2", ",", "a2", ",", 
       RowBox[{"p", "+", "a4"}], ",", "0", ",", "0"}], "}"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"A\[Gamma]", "=", 
     RowBox[{"{", 
      RowBox[{"0", ",", "0", ",", 
       RowBox[{"-", "b4"}], ",", "p", ",", "0"}], "}"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"A\[Psi]", "=", 
     RowBox[{"{", 
      RowBox[{"0", ",", "0", ",", "0", ",", "0", ",", "p"}], "}"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{"MatrixForm", "[", 
    RowBox[{"Abok", "=", 
     RowBox[{"Transpose", "[", 
      RowBox[{"{", 
       RowBox[{
       "A\[Omega]x", ",", "A\[Omega]y", ",", "A\[Beta]", ",", "A\[Gamma]", 
        ",", "A\[Psi]"}], "}"}], "]"}]}], "]"}]}]}]], "Input"],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
   "\:0412\:0435\:043a\:0442\:043e\:0440", " ", "U", " ", 
    "\:043f\:0440\:043e\:0434\:043e\:043b\:044c\:043d\:043e\:0433\:043e", " ",
     "\:0434\:0432\:0438\:0436\:0435\:043d\:0438\:044f"}], "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{"Ubok", "=", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"\[Delta]\:043d", "[", "t", "]"}], ",", 
     RowBox[{"\[Delta]\:044d", "[", "t", "]"}]}], "}"}]}]}]], "Input"],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
   "\:041c\:0430\:0442\:0440\:0438\:0446\:0430", " ", "B", " ", 
    "\:043f\:0440\:043e\:0434\:043e\:043b\:044c\:043d\:043e\:0433\:043e", " ",
     "\:0434\:0432\:0438\:0436\:0435\:043d\:0438\:044f"}], "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"B\[Delta]\:043d", "=", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"-", "a5"}], ",", 
       RowBox[{"-", "a3"}], ",", 
       RowBox[{"-", "a7"}], ",", "0", ",", "0"}], "}"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"B\[Delta]\:044d", "=", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"-", "b3"}], ",", 
       RowBox[{"-", "b5"}], ",", "0", ",", "0", ",", "0"}], "}"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{"MatrixForm", "[", 
    RowBox[{"Bbok", "=", 
     RowBox[{"Transpose", "[", 
      RowBox[{"{", 
       RowBox[{"B\[Delta]\:043d", ",", "B\[Delta]\:044d"}], "}"}], "]"}]}], 
    "]"}]}]}]], "Input"],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
   "\:0423\:0440\:0430\:0432\:043d\:0435\:043d\:0438\:044f", " ", 
    "\:043f\:0440\:043e\:0434\:043e\:043b\:044c\:043d\:043e\:0433\:043e", " ",
     "\:0434\:0432\:0438\:0436\:0435\:043d\:0438\:044f"}], "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{"MatrixForm", "[", 
   RowBox[{"prod", "=", 
    RowBox[{"Table", "[", " ", "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{
       RowBox[{"Sum", "[", 
        RowBox[{
         RowBox[{
          RowBox[{"Abok", "[", 
           RowBox[{"[", 
            RowBox[{"i", ",", "j"}], "]"}], "]"}], "*", 
          RowBox[{"Xbok", "[", 
           RowBox[{"[", "j", "]"}], "]"}]}], ",", 
         RowBox[{"{", 
          RowBox[{"j", ",", "1", ",", "5"}], "}"}]}], "]"}], "\[Equal]", 
       RowBox[{"Sum", "[", 
        RowBox[{
         RowBox[{
          RowBox[{"Bbok", "[", 
           RowBox[{"[", 
            RowBox[{"i", ",", "j"}], "]"}], "]"}], "*", 
          RowBox[{"Ubok", "[", 
           RowBox[{"[", "j", "]"}], "]"}]}], ",", 
         RowBox[{"{", 
          RowBox[{"j", ",", "1", ",", "2"}], "}"}]}], "]"}]}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", "1", ",", "5"}], "}"}]}], " ", "]"}]}], 
   "]"}]}]], "Input"],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
   "\:0417\:0430\:0434\:0430\:0434\:0438\:043c\:0441\:044f", " ", 
    "\:043a\:043e\:044d\:0444\:0444\:0438\:0446\:0438\:0435\:043d\:0442\:0430\
\:043c\:0438"}], "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"c1", "=", "0.587"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"c2", "=", "45.8"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"c3", "=", "13.6"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"c4", "=", "1.21"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"c5", "=", "0.204"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"c6", "=", "5.721"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"a1", "=", "1.64"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"a2", "=", "40"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"a3", "=", "7.56"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"a4", "=", "0.62"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"a5", "=", "7.36"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"a6", "=", "0.985"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"a7", "=", "0.0278"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{
     RowBox[{"a8", "=", "1"}], ";"}], "*)"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"b1", "=", "7.235"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"b2", "=", "133"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"b3", "=", "22.9"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"b4", "=", "0.0191"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"b5", "=", 
     RowBox[{"-", "1.05"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"b6", "=", "0.131"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"b7", "=", "0.004"}], ";"}]}]}]], "Input"],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
   "\:041e\:0446\:0435\:043d\:0438\:043c", " ", 
    "\:0432\:043e\:0437\:043c\:043e\:0436\:043d\:043e\:0441\:0442\:044c", " ",
     "\:0438\:0437\:043e\:043b\:0438\:0440\:043e\:0432\:0430\:043d\:043d\:043e\
\:0433\:043e", " ", 
    "\:0438\:0441\:0441\:043b\:0435\:0434\:043e\:0432\:0430\:043d\:0438\:044f\
", " ", "\:0434\:0432\:0438\:0436\:0435\:043d\:0438\:044f", " ", 
    "\"\<\:0440\:044b\:0441\:043a\:0430\:043d\:0438\:0435-\:0441\:043a\:043e\
\:043b\:044c\:0436\:0435\:043d\:0438\:0435\>\"", " ", "\:0438", " ", 
    "\:0434\:0432\:0438\:0436\:0435\:043d\:0438\:044f", " ", "\:043f\:043e", 
    " ", "\:043a\:0440\:0435\:043d\:0443"}], "*)"}], "\[IndentingNewLine]", 
  RowBox[{"If", "[", " ", 
   RowBox[{
    RowBox[{
     RowBox[{"b1", " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"a1", " ", "a4"}], "+", "a2"}], ")"}]}], "\[GreaterEqual]", 
     RowBox[{"9", " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"a1", " ", "b2", " ", "b7"}], "-", 
        RowBox[{"b2", " ", "b6"}], "-", 
        RowBox[{"a6", " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"a4", " ", "b6"}], "+", 
           RowBox[{"a2", " ", "b1"}]}], ")"}]}]}], ")"}]}]}], ",", 
    RowBox[{"Print", "[", "\"\<\:041c\:043e\:0436\:043d\:043e\>\"", "]"}], 
    ",", 
    RowBox[{
    "Print", "[", "\"\<\:041d\:0435\:043b\:044c\:0437\:044f\>\"", "]"}]}], 
   "]"}]}]], "Input"],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
   "\:041a\:043e\:044d\:0444\:0438\:0446\:0438\:0435\:043d\:0442", " ", 
    "\:0434\:0435\:043c\:043f\:0444\:0438\:0440\:043e\:0432\:0430\:043d\:0438\
\:044f"}], " ", "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"\[Xi]\[Beta]", "=", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{"a1", "+", "a4"}], ")"}], "/", 
     RowBox[{"(", 
      RowBox[{"2", "*", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{"a2", "+", 
          RowBox[{"a1", " ", "a4"}]}], ")"}], "^", "0.5"}]}], ")"}]}]}], 
   "\[IndentingNewLine]", 
   RowBox[{"If", "[", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{
       RowBox[{"\[Xi]\[Beta]", "<", "0.2"}], " ", "||", " ", 
       RowBox[{"\[Xi]\[Beta]", ">", "0.4"}]}], ")"}], ",", 
     RowBox[{
     "Print", "[", 
      "\"\<\:041d\:0443\:0436\:0435\:043d \
\:0434\:0435\:043c\:043f\:0444\:0435\:0440 \:0440\:044b\:0441\:043a\:0430\
\:043d\:0438\:044f\>\"", "]"}], ",", 
     RowBox[{
     "Print", "[", 
      "\"\<\:041d\:0435 \:043d\:0443\:0436\:0435\:043d \:0434\:0435\:043c\
\:043f\:0444\:0435\:0440 \
\:0440\:044b\:0441\:043a\:0430\:043d\:0438\:044f\>\"", "]"}]}], 
    "]"}]}]}]], "Input"],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
   "\:0427\:0430\:0441\:0442\:043e\:0442\:0430", " ", 
    "\:0441\:043e\:0431\:0441\:0442\:0432\:0435\:043d\:043d\:044b\:0445", " ",
     "\:043a\:043e\:043b\:0435\:0431\:0430\:043d\:0438\:0439", " ", 
    "\:0441\:0432\:043e\:0431\:043e\:0434\:043d\:043e\:0433\:043e", " ", 
    "\:0434\:0432\:0438\:0436\:0435\:043d\:0438\:044f", " ", "\:043f\:043e", 
    " ", "\:0443\:0433\:043b\:0443", " ", 
    "\:0440\:044b\:0441\:043a\:0430\:043d\:0438\:044f"}], "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"T\[Beta]", "=", 
    RowBox[{"1", "/", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"a2", "+", 
        RowBox[{"a1", " ", "a4"}]}], ")"}], "^", "0.5"}]}]}], 
   "\[IndentingNewLine]", 
   RowBox[{"\[CapitalOmega]\[Beta]", "=", 
    RowBox[{"1", "/", "T\[Beta]"}]}]}]}]], "Input"],

Cell[TextData[StyleBox["\:0418\:0437\:043e\:043b\:0438\:0440\:043e\:0432\:0430\
\:043d\:043d\:043e\:0435 \:0434\:0432\:0438\:0436\:0435\:043d\:0438\:0435 \
\:0441\:0430\:043c\:043e\:043b\:0435\:0442\:0430 \:0432 \:0443\:0433\:043b\
\:0430\:0445 \"\:0440\:044b\:0441\:043a\:0430\:043d\:0438\:0435-\:0441\:043a\
\:043e\:043b\:044c\:0436\:0435\:043d\:0438\:0435\"",
 FontFamily->"Times New Roman",
 FontSize->48,
 FontColor->RGBColor[0., 0., 1.]]], "Text"],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
   "\:0412\:0435\:043a\:0442\:043e\:0440", " ", "X", " ", 
    "\:0431\:043e\:043a\:043e\:0432\:043e\:0433\:043e", " ", 
    "\:0434\:0432\:0438\:0436\:0435\:043d\:0438\:044f"}], "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{"Xbok1", "=", 
   RowBox[{"{", " ", 
    RowBox[{
     RowBox[{"\[Omega]y", "[", "t", "]"}], ",", 
     RowBox[{"\[Beta]", "[", "t", "]"}]}], "}"}]}]}]], "Input"],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
   "\:041c\:0430\:0442\:0440\:0438\:0446\:0430", " ", "\:0410", " ", 
    "\:0431\:043e\:043a\:043e\:0432\:043e\:0433\:043e", " ", 
    "\:0434\:0432\:0438\:0436\:0435\:043d\:0438\:044f"}], "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"A\[Omega]y", "=", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"p", "+", "a1"}], ",", 
       RowBox[{"-", "1"}]}], "}"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"A\[Beta]", "=", 
     RowBox[{"{", 
      RowBox[{"a2", ",", 
       RowBox[{"p", "+", "a4"}]}], "}"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{
     RowBox[{"A\[Omega]y", "=", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"p", "+", "a1"}], ",", "a2"}], "}"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"A\[Beta]", "=", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"-", "a8"}], ",", 
        RowBox[{"p", "+", "a4"}]}], "}"}]}], ";"}], "*)"}], 
   "\[IndentingNewLine]", 
   RowBox[{"MatrixForm", "[", 
    RowBox[{"Abok1", "=", 
     RowBox[{"Transpose", "[", 
      RowBox[{"{", 
       RowBox[{"A\[Omega]y", ",", "A\[Beta]"}], "}"}], "]"}]}], 
    "]"}]}]}]], "Input"],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
   "\:0412\:0435\:043a\:0442\:043e\:0440", " ", "U", " ", 
    "\:043f\:0440\:043e\:0434\:043e\:043b\:044c\:043d\:043e\:0433\:043e", " ",
     "\:0434\:0432\:0438\:0436\:0435\:043d\:0438\:044f"}], "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{"Ubok1", "=", 
   RowBox[{"{", 
    RowBox[{"\[Delta]\:043d", "[", "t", "]"}], "}"}]}]}]], "Input"],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
   "\:041c\:0430\:0442\:0440\:0438\:0446\:0430", " ", "B", " ", 
    "\:043f\:0440\:043e\:0434\:043e\:043b\:044c\:043d\:043e\:0433\:043e", " ",
     "\:0434\:0432\:0438\:0436\:0435\:043d\:0438\:044f"}], "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"B\[Delta]\:043d", "=", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"-", "a3"}], ",", 
       RowBox[{"-", "a7"}]}], "}"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{"MatrixForm", "[", 
    RowBox[{"Bbok1", "=", 
     RowBox[{"{", "B\[Delta]\:043d", "}"}]}], "]"}]}]}]], "Input"],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
   "\:0423\:0440\:0430\:0432\:043d\:0435\:043d\:0438\:044f", " ", 
    "\:043f\:0440\:043e\:0434\:043e\:043b\:044c\:043d\:043e\:0433\:043e", " ",
     "\:0434\:0432\:0438\:0436\:0435\:043d\:0438\:044f"}], "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{"MatrixForm", "[", 
   RowBox[{"prod", "=", 
    RowBox[{"Table", "[", " ", 
     RowBox[{
      RowBox[{
       RowBox[{"Sum", "[", 
        RowBox[{
         RowBox[{
          RowBox[{"Abok1", "[", 
           RowBox[{"[", 
            RowBox[{"i", ",", "j"}], "]"}], "]"}], "*", 
          RowBox[{"Xbok1", "[", 
           RowBox[{"[", "j", "]"}], "]"}]}], ",", 
         RowBox[{"{", 
          RowBox[{"j", ",", "1", ",", "2"}], "}"}]}], "]"}], "\[Equal]", 
       RowBox[{
        RowBox[{"Bbok1", "[", 
         RowBox[{"[", 
          RowBox[{"1", ",", "i"}], "]"}], "]"}], "*", 
        RowBox[{"Ubok1", "[", 
         RowBox[{"[", "1", "]"}], "]"}]}]}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", "1", ",", "2"}], "}"}]}], " ", "]"}]}], 
   "]"}]}]], "Input"],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
   "\:0417\:0430\:043f\:0438\:0448\:0435\:043c", " ", 
    "\:043f\:0435\:0440\:0435\:0434\:0430\:0442\:043e\:0447\:043d\:0443\:044e\
", " ", "\:0444\:0443\:043d\:043a\:0446\:0438\:044e", " ", 
    RowBox[{"Wy", "/", "\[Delta]\:043d"}], " ", 
    "\:0437\:0430\:043c\:043a\:043d\:0443\:0442\:043e\:0439", " ", 
    "\:0441\:0438\:0441\:0442\:0435\:043c\:044b"}], "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{"W0z", "=", 
   RowBox[{"Simplify", "[", 
    RowBox[{
     RowBox[{"Det", "[", 
      RowBox[{"{", 
       RowBox[{"B\[Delta]\:043d", ",", "A\[Beta]"}], "}"}], "]"}], "/", 
     RowBox[{"Det", "[", 
      RowBox[{"{", 
       RowBox[{"A\[Omega]y", ",", "A\[Beta]"}], "}"}], "]"}]}], 
    "]"}]}]}]], "Input"],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
   "\:041f\:043e\:0441\:0442\:0440\:043e\:0438\:043c", " ", 
    "\:043f\:0435\:0440\:0435\:0434\:0430\:0442\:043e\:0447\:043d\:0443\:044e\
", " ", "\:0444\:0443\:043d\:043a\:0446\:0438\:044e", " ", 
    RowBox[{"Wy", "/", "\[Delta]\:043d"}], " ", 
    "\:0437\:0430\:043c\:043a\:043d\:0443\:0442\:043e\:0439", " ", 
    "\:0441\:0438\:0441\:0442\:0435\:043c\:044b"}], "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{"Step", "[", 
   RowBox[{"W0z", ",", "6"}], "]"}]}]], "Input"],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
   "\:0417\:0430\:043f\:0438\:0448\:0435\:043c", " ", 
    "\:043f\:0435\:0440\:0435\:0434\:0430\:0442\:043e\:0447\:043d\:0443\:044e\
", " ", "\:0444\:0443\:043d\:043a\:0446\:0438\:044e", " ", 
    RowBox[{"Wy", "/", "\[Delta]\:043d"}], " ", 
    "\:0440\:0430\:0437\:043e\:043c\:043a\:043d\:0443\:0442\:043e\:0439", " ",
     "\:0441\:0438\:0441\:0442\:0435\:043c\:044b"}], "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{"W0r", "=", 
   RowBox[{"Simplify", "[", 
    FractionBox["W0z", 
     RowBox[{"1", "-", "W0z"}]], "]"}]}]}]], "Input"],

Cell[BoxData[
 RowBox[{"\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{
   "\:041f\:043e\:0441\:0442\:0440\:043e\:0438\:043c", " ", 
    "\:041b\:0410\:0427\:0425", " ", 
    RowBox[{"Wy", "/", "\[Delta]\:043d"}], " ", 
    "\:0440\:0430\:0437\:043e\:043c\:043a\:043d\:0443\:0442\:043e\:0439", " ",
     "\:0441\:0438\:0441\:0442\:0435\:043c\:044b"}], "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{"Bode", "[", 
   RowBox[{"W0r", ",", 
    RowBox[{"-", "1"}], ",", "3"}], "]"}]}]], "Input"],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
   "\:0420\:0430\:0441\:0441\:0447\:0438\:0442\:0430\:0435\:043c", " ", 
    "\:043f\:0430\:0440\:0430\:043c\:0435\:0442\:0440\:044b", " ", 
    "\:0434\:0435\:043c\:043f\:0444\:0435\:0440\:0430", " ", 
    "\:0440\:044b\:0441\:043a\:0430\:043d\:0438\:044f"}], "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"\[Mu]", "=", 
    FractionBox[
     RowBox[{
      RowBox[{"2", "*", " ", 
       SqrtBox[
        RowBox[{"a2", "+", " ", 
         RowBox[{"a1", " ", "a4"}]}]]}], "-", 
      RowBox[{"(", 
       RowBox[{"a1", "+", "a4"}], ")"}]}], "a3"]}], "\[IndentingNewLine]", 
   "\[IndentingNewLine]", 
   RowBox[{"Tn", "=", 
    FractionBox["a3", 
     RowBox[{
      RowBox[{"a3", " ", "a4"}], "-", 
      RowBox[{"a2", " ", "a7"}]}]]}]}]}]], "Input"],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
   "\:041e\:0431\:0440\:0430\:0442\:043d\:0430\:044f", " ", 
    "\:0441\:0432\:044f\:0437\:044c"}], "*)"}], "\[IndentingNewLine]", 
  RowBox[{"Wu", "=", 
   RowBox[{"Simplify", "[", 
    RowBox[{"\[Mu]", "*", 
     FractionBox[
      RowBox[{"Tn", "*", "p"}], 
      RowBox[{
       RowBox[{"Tn", "*", "p"}], "+", "1"}]]}], "]"}]}]}]], "Input"],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
   "\:0417\:0430\:043f\:0438\:0448\:0435\:043c", " ", 
    "\:043f\:0435\:0440\:0435\:0434\:0430\:0442\:043e\:0447\:043d\:0443\:044e\
", " ", "\:0441\:043a\:043e\:0440\:0435\:043a\:0442\:0438\:0440\:043e\:0432\
\:0430\:043d\:043d\:0443\:044e", " ", 
    "\:0444\:0443\:043d\:043a\:0446\:0438\:044e", " ", 
    "\:0437\:0430\:043c\:043a\:043d\:0443\:0442\:043e\:0439", " ", 
    "\:0441\:0438\:0441\:0442\:0435\:043c\:044b", " ", 
    RowBox[{"Wy", "/", "\[Delta]\:043d"}]}], "*)"}], "\[IndentingNewLine]", 
  RowBox[{"W1z", "=", 
   RowBox[{"Simplify", "[", 
    FractionBox["W0z", 
     RowBox[{"1", "-", 
      RowBox[{"W0z", "*", "Wu"}]}]], "]"}]}]}]], "Input"],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
   "\:041f\:043e\:0441\:0442\:0440\:043e\:0438\:043c", " ", 
    "\:043f\:0435\:0440\:0435\:0434\:0430\:0442\:043e\:0447\:043d\:0443\:044e\
", " ", "\:0437\:0430\:043c\:043a\:043d\:0443\:0442\:043e\:0439", " ", 
    "\:0441\:043a\:043e\:0440\:0435\:043a\:0442\:0438\:0440\:043e\:0432\:0430\
\:043d\:043d\:043e\:0439", " ", "\:0444\:0443\:043d\:043a\:0446\:0438\:044e", 
    " ", 
    RowBox[{"Wy", "/", "\[Delta]\:043d"}]}], "*)"}], "\[IndentingNewLine]", 
  RowBox[{"Step", "[", 
   RowBox[{"W1z", ",", "6"}], "]"}]}]], "Input"],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
   "\:0417\:0430\:043f\:0438\:0448\:0435\:043c", " ", 
    "\:043f\:0435\:0440\:0435\:0434\:0430\:0442\:043e\:0447\:043d\:0443\:044e\
", " ", "\:0444\:0443\:043d\:043a\:0446\:0438\:044e", " ", 
    RowBox[{"Wy", "/", "\[Delta]\:043d"}], " ", 
    "\:0440\:0430\:0437\:043e\:043c\:043a\:043d\:0443\:0442\:043e\:0439", " ",
     "\:0441\:043a\:043e\:0440\:0435\:043a\:0442\:0438\:0440\:043e\:0432\:0430\
\:043d\:043d\:043e\:0439", " ", 
    "\:0441\:0438\:0441\:0442\:0435\:043c\:044b"}], "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{"W1r", "=", 
   RowBox[{"Simplify", "[", 
    FractionBox["W1z", 
     RowBox[{"1", "-", "W1z"}]], "]"}]}]}]], "Input"],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
   "\:041f\:043e\:0441\:0442\:0440\:043e\:0438\:043c", " ", 
    "\:041b\:0410\:0427\:0425", " ", 
    RowBox[{"Wy", "/", "\[Delta]\:043d"}], " ", 
    "\:0440\:0430\:0437\:043e\:043c\:043a\:043d\:0443\:0442\:043e\:0439", " ",
     "\:0441\:043a\:043e\:0440\:0435\:043a\:0442\:0438\:0440\:043e\:0432\:0430\
\:043d\:043d\:043e\:0439", " ", 
    "\:0441\:0438\:0441\:0442\:0435\:043c\:044b"}], "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{"Bode", "[", 
   RowBox[{"W1r", ",", 
    RowBox[{"-", "1"}], ",", "3"}], "]"}]}]], "Input"]
},
WindowSize->{907, 701},
WindowMargins->{{2, Automatic}, {Automatic, 0}},
CellContext->Notebook,
TrackCellChangeTimes->False,
FrontEndVersion->"7.0 for Microsoft Windows (32-bit) (February 18, 2009)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[545, 20, 362, 6, 127, "Text"],
Cell[910, 28, 73, 1, 31, "Input"],
Cell[986, 31, 1968, 52, 152, "Input"],
Cell[2957, 85, 514, 14, 31, "Input"],
Cell[3474, 101, 228, 4, 69, "Text"],
Cell[3705, 107, 567, 14, 52, "Input"],
Cell[4275, 123, 1380, 40, 132, "Input"],
Cell[5658, 165, 394, 10, 52, "Input"],
Cell[6055, 177, 443, 12, 52, "Input"],
Cell[6501, 191, 1078, 30, 72, "Input"],
Cell[7582, 223, 210, 4, 69, "Text"],
Cell[7795, 229, 581, 15, 52, "Input"],
Cell[8379, 246, 1586, 46, 152, "Input"],
Cell[9968, 294, 463, 12, 52, "Input"],
Cell[10434, 308, 978, 29, 92, "Input"],
Cell[11415, 339, 1248, 35, 72, "Input"],
Cell[12666, 376, 1898, 51, 452, "Input"],
Cell[14567, 429, 1451, 35, 52, "Input"],
Cell[16021, 466, 1195, 35, 72, "Input"],
Cell[17219, 503, 852, 21, 72, "Input"],
Cell[18074, 526, 453, 7, 127, "Text"],
Cell[18530, 535, 439, 12, 52, "Input"],
Cell[18972, 549, 1217, 38, 132, "Input"],
Cell[20192, 589, 394, 10, 52, "Input"],
Cell[20589, 601, 613, 17, 72, "Input"],
Cell[21205, 620, 1082, 31, 52, "Input"],
Cell[22290, 653, 770, 20, 52, "Input"],
Cell[23063, 675, 529, 12, 52, "Input"],
Cell[23595, 689, 592, 14, 67, "Input"],
Cell[24190, 705, 493, 12, 72, "Input"],
Cell[24686, 719, 816, 24, 131, "Input"],
Cell[25505, 745, 398, 12, 69, "Input"],
Cell[25906, 759, 717, 16, 67, "Input"],
Cell[26626, 777, 582, 12, 52, "Input"],
Cell[27211, 791, 704, 16, 67, "Input"],
Cell[27918, 809, 582, 14, 52, "Input"]
}
]
*)

(* End of internal cache information *)